Enhanced Kerr nonlinearity in a quantized four-level graphene nanostructure
Solookinejad Ghahraman, Panahi M, Ahmadi E, Hossein Asadpour Seyyed†,
Department of Physics, Marvdasht Branch, Islamic Azad University, Marvdasht, Iran

 

† Corresponding author. E-mail: S.Hosein.Asadpour@gmail.com

Abstract
Abstract

In this paper, a new model is proposed for manipulating the Kerr nonlinearity of right-hand circular probe light in a monolayer of graphene nanostructure. By using the density matrix equations and quantum optical approach, the third-order susceptibility of probe light is explored numerically. It is realized that the enhanced Kerr nonlinearity with zero linear absorption can be provided by selecting the appropriate quantities of controllable parameters, such as Rabi frequency and elliptical parameter of elliptical polarized coupling field. Our results may be useful applications in future all-optical system devices in nanostructures.

1. Introduction

Electromagnetically induced transparency (EIT),[1] in atomic systems and semiconductor quantum wells and quantum dot structures has led to many unforeseen phenomena such as lasing without inversion,[2] four-wave mixing,[3,4] and so on.[515] EIT exploits destructive quantum interference to make an otherwise opaque medium almost transparent. This feature of EIT medium (transparent to the applied fields) makes it suitable for application in nonlinear optical science. In nonlinear optics, it is favorable to have a strong nonlinear medium in low light intensities. Therefore, Kerr nonlinearity of such medium plays essential roles for studying some other phenomena which require strong nonlinear medium such as quantum nondemolition measurements,[16] quantum bit regeneration,[17] quantum state teleportation,[18] and the generation of the optical solitons.[19] The Kerr nonlinearity of the medium corresponds to the real part of third-order susceptibility of weak probe light which propagates through it. The features of Kerr nonlinearity of weak probe light in multi-level atomic systems have been reported by many research groups.[2029] It has been shown that the Kerr nonlinearity can be manipulated by many different mechanisms. For instance, Niu et al. proposed a model for giant Kerr-nonlinearity-based spontaneously generated coherence (SGC) in a usual three-level atomic system.[26] In another study done by Niu et al., the role of the double dark resonances on Kerr nonlinearity has been discussed.[27] Asadpour et al.[28] proposed a model for studying the enhanced Kerr nonlinearity via quantum interference from spontaneous emission. The properties of Kerr nonlinearity in a semiconductor quantum well and quantum dot-based different mechanisms have also been investigated.[3035] For instance, phase control of Kerr nonlinearity in a quantum well nanostructure has been suggested by Asadpour et al.[30] The role of biexciton coherence on Kerr nonlinearity of probe light propagated in a multiple quantum well nanostructure has also been reported.[33] Moreover, it has been found that exciton spin relaxation can lead to manipulating of Kerr nonlinearity in a multiple quantum well nanostructure.[34] The giant Kerr nonlinearity in a crystal of molecular magnet has also been discussed.[35]

On the other hand, the optical and electronic features of graphene nanostructure due to its important usages in condense matter and optical physics have been investigated recently.[3639] The properties of electrons near the Dirac points make it possible for interacting by electromagnetic fields.[4042] The nonlinear optical features of graphene nanostructure in the presence of strong magnetic field has been discussed very recently in four-wave mixing,[43,44] generation of polarized-entangle photon, nonlinear frequency conversion of THz surface plasmons based on the nonlinear optical interaction,[45,47] optical bistability,[48,49] and slow and fast light propagating.[50] To the best of our knowledge, the enhanced Kerr nonlinearity of graphene with zero linear and nonlinear absorption has not been reported theoretically or experimentally. In this paper, we suggest a model based on the quantized four-level graphene nanostructure for enhancing the Kerr nonlinearity by using a probe and an elliptical polarized control field. We show that by selecting appropriate values for Rabi frequency and elliptical parameter of control field, the Kerr nonlinearity can be enhanced at certain values of probe field detuning.

2. Model and equations

A single layer of graphene system in the presence of strong magnetic field is presented in a four-level scheme (Fig. 1). Due to the presence of external magnetic field, the unbroken energy bands near Dirac points split into individual Landau levels. It is well known that the selection rules in graphene obey from Δ|n| = ±1, where n is the energy quantum number. It should be noted that the optical transitions between Landau levels in the range of 0.01 T–10 T can be located in the infrared and terahertz regions due to ħωc ≃ 36 B1/2·T1/2·meV. The mentioned graphene system interacts with a probe and an elliptical polarized control field. The elliptically polarized control field can be regarded as a combination of the right- and left-circularly polarized components which can be obtained by using a quarter wave plate (QWP).[51] After passing a vertical polarized coupling field with amplitude E0 through QWP, the elliptical polarized coupling field can be obtained with E0 = E+ σ+ + Eσ, where , and . Parameter θ is the ellipticity parameter. The unit vectors of right- and left-circularly polarized basis are denoted by σ+ and σ, respectively. The probe light with right-hand circular polarization, frequency ωp, and Rabi frequency Ωp drives transition |1⟩ ↔ |4⟩ and elliptical polarized coupling field with frequency ω0 and Rabi frequency Ωc act on transitions |4⟩ ↔ |3⟩ and |3⟩ ↔ |2⟩, respectively. The Rabi frequencies for transitions |4⟩ ↔ |3⟩ and |3⟩ ↔ |2⟩ are equal to Ω3 = Ωc(cos θ − sin θ)e−iθ and Ω3 = Ωc (cos θ + sin θ)eiθ, respectively. It is assumed that μ43 = μ32 = μ and . The corresponding optical transition between |m⟩ and |n⟩ denoted by μ = ⟨m|μ|n⟩ = em|r|n⟩ = iħ/εmεnm|vfσ|n⟩, where σ = (σx,σy) is the Pauli matrix vector and vF = 3γ0/2ħa ≃ 106 m/s is Fermi velocity (γ0 ≃ 2.8 eV and a = 1.42 Å are the nearest neighbor hopping energy and C–C spacing).[52] The corresponding Hamiltonian for a monolayer of graphene system is governed by

where denotes the generalized momentum operator, is the electron momentum operator, e is the electron charge, and A is the vector potential, which is equal to (0, Bx) for a static magnetic field. For considering the incident optical field on the Hamiltonian of the system, one needs to add the vector potential of the optical field (Aopt = icE/ω and E = Ep + E0) to the vector potential of the magnetic field in the generalized momentum operator . After this interchange, the interaction Hamiltonian can be written in the following form:

A standard dynamic equation for the density matrix of Dirac electrons in graphene by using the Liouville’s equation

can be obtained. Here, indicates incoherent relaxation which may originate from disorder, interaction with phonons, and carrier–carrier interactions. The density matrix elements of the system are governed by

where Δp = (εn = 1εn = − 2)/ħωp, Δ2 = (εn = 0εn = −1)/ħω0, and Δ3 = (εn = 1εn = 0)/ħω0 represent the corresponding frequency detunings, and is the energy of the Landau level for electrons near the Dirac point, with , and implies the magnetic length. For achieving the third-order of susceptibility, we assume that the probe field is weaker than the elliptical coupling field. In this case density matrix elements are expanding as . The zeroth order solution of ρ11 will be identical, i.e., , and other elements are set to be zero. The first- and third-order susceptibilities χ(1) and χ(3) of the medium can be determined by coherence terms and , respectively, as

where N and εr are the sheet electron density of graphene and the substrate dielectric constant, respectively.

Fig. 1. (a) The LLs near the K point superimposed on the electronic energy dispersion without a magnetic field E = ±vF|p|. The magnetic field condenses the original states in the Dirac cone into discrete energies. The LLs in graphene are unequally spaced: . (b) Energy level diagram and optical transitions in a monolayer graphene interacting with probe and elliptical polarized coupling field. Graphene monolayer is a one-atom-thick monolayer of carbon atoms arranged in a hexagonal lattice, which we will treat as a perfect two-dimensional (2D) crystal structure in the xy plane.
3. Results and discussion

In the following section, we will analyze the third-order susceptibility through numerical simulations in detail. According to the numerical estimate based on Refs. [43]–[50], we can take a reasonable value for the decay rate and all parameters are scaled by decay rate γ. This value depends on the sample quality and the substrate used in experiments.[51,52] The electron concentration is N ≃ 5 × 1012 cm−2 and the substrate dielectric constant is εr ≃ 4.5.[5355] The linear and nonlinear features of the susceptibility for parameter θ = 0 are explored in Fig. 2. By choosing θ = 0, the Rabi-frequencies of elliptical polarized field will become Ω+ = Ω = Ωc. In this case, we will expect that due to the double dark resonances, the medium is accompanied with strong linear and nonlinear absorption. From Fig. 2, one can find that the Kerr nonlinearity of the medium (dashed line of Fig. 2(b)) corresponds to strong linear and nonlinear absorption (solid line of Figs. 2(a) and 2(b)). These properties are not suitable for using in nonlinear optics where enhanced Kerr nonlinearity with zero linear absorption is needed. For overcoming this deficiency, we examine the effect of elliptical parameter θ on the behaviors of linear and nonlinear susceptibilities. The properties of linear and nonlinear susceptibilities for θ = π/6 are presented in Fig. 3. The linear properties are shown in Fig. 3(a) and nonlinear features are presented in Fig. 3(b). It can be easily seen that the linear absorption (solid line of Fig. 3(a)) of weak probe light at Δp = ±0.41γ is zero. This means that the medium for these quantities of detuning is transparent for the probe light. The slope of dispersion curve in this situation is positive corresponding to slow light propagation. Until now, these features of linear susceptibility are favorable for applications in nonlinear optics. The nonlinear behaviors of susceptibility for θ = π/6 are presented in Fig. 3(b). It can be realized that the nonlinear absorption converts to nonlinear amplification in two detunings of probe field. In this situation, the Kerr nonlinearity is enhanced and reaches to maximum quantity at Δp = ±0.41γ. These properties of the medium due to zero linear absorption and enhanced Kerr nonlinearity are very suitable for application in nonlinear optical technologies, where the enhanced Kerr nonlinearity at low light level is needed. Compared with the single dark resonance case, the interacting double dark resonances cause an enhancement of the Kerr nonlinearity with vanishing linear absorption. In the dressed-state picture, this case corresponds to the profound interference situation because the eigenvalues of two of the dressed states intersect.[56] For further discussion about the role of elliptical parameter θ on linear absorption and Kerr nonlinearity of the medium, we display the features of linear absorption and Kerr nonlinearity of medium versus θ in Fig. 4. It is realized that the enhanced Kerr nonlinearity is accompanied by vanishing absorption at different values of parameter θ. These results have good agreement with our previous results. At the final step, the properties of linear absorption (solid line) and Kerr nonlinearity (dashed line) versus Rabi frequency Ωc for θ = 0 (Fig. 4(a)), and θ = π/6 (Fig. 4(b)) are plotted in Fig. 5. By careful attention to the obtained results, one can find that when we choose θ = 0, by enhancing the Rabi frequency Ωc, the Kerr nonlinearity increases and linear absorption reduces, however, the linear absorption never reaches zero. In addition, when we select θ = π/6, the Kerr nonlinearity is enhanced and simultaneously linear absorption is vanished when Ωc reaches about 0.4γ. This is very interesting for us to reach enhanced Kerr nonlinearity with vanishing absorption. As a result, we find that the enhanced Kerr nonlinearity with vanishing absorption is caused by simultaneous effects of elliptical parameter and Rabi frequency Ωc. In fact, the elliptical parameter θ leads to modifying intensities of the Rabi frequencies Ω2 and Ω3, respectively. This means that the enhanced Kerr nonlinearity with zero linear absorption can be tuned by elliptical parameter.

Fig. 2. (a) Linear susceptibility χ(1) versus probe light detuning. Solid line corresponds to linear absorption and dashed line corresponds to dispersion. (b) Nonlinear susceptibility χ(3) versus probe light detuning. Solid line corresponds to nonlinear absorption and dashed line corresponds to Kerr nonlinearity. The used parameters are Ωp = 0.01γ, Ωc = 0.4γ, Δ2 = Δ3 = 0, and θ = 0.
Fig. 3. (a) Linear susceptibility χ(1) versus probe light detuning. Solid line corresponds to linear absorption and dashed line corresponds to dispersion. (b) Nonlinear susceptibility χ(3) versus probe light detuning. Solid line corresponds to nonlinear absorption and dashed line corresponds to Kerr nonlinearity. Parameter θ = π/6 and others are the same as Fig. 2.
Fig. 4. Linear absorption (solid line) and Kerr nonlinearity (dashed line) versus elliptical parameter θ. The selected parameters are the same as Fig. 2.
Fig. 5. Linear absorption (solid line) and Kerr nonlinearity (dashed line) versus different quantities of Rabi frequency Ωc for (a) θ = 0 and (b) θ = π/6. The other parameters are the same as Fig. 2.

Similar to the laser-driven atomic medium or ensemble consisting of many single atoms,[57] here we only consider a graphene monolayer in the graphene ensemble interacting with optical fields for the sake of simplification of the calculation. It should be emphasized that the optical conductivity of a graphene layer is proportional to its effective thickness or surface density of atoms.[58,59] This is in spite of the fact that the effective thickness of a graphene layer is around one angstrom and the vacuum wavelength of the optical waves under consideration is around one micron. However, through our calculation and analysis, it is possible to observe such spatial solitons by use of the graphene ensemble. In view of rapid advances in graphene material, we believe that the formation of matched infrared soliton pairs will be accessible in experiments in the near future.

4. Conclusion

In summary, the linear and nonlinear susceptibilities of the monolayer graphene system driven by a probe and elliptical polarized coupling fields are explored by using the density matrix formalism. Our numerical results show that by tuning the elliptical parameter of elliptical polarized coupling field, the Rabi frequencies of coupling fields can be modified. Therefore, the enhanced Kerr nonlinearity with vanishing of linear absorption can be provided in two detunings of probe field.

Reference
1Wu YYang X 2005 Phys. Rev. 71 053806
2Zhao YHuang DWu C 1996 J. Opt. Soc. Am. B: Opt. Phys. 13 1614
3Wu YDeng L 2004 Phys. Rev. Lett. 93 143904
4Wu YYang X 2004 Phys. Rev. 70 053818
5Wang ZJiang J 2010 Phys. Lett. 374 4853
6Wu YYang X 2005 Phys. Rev. 71 053806
7Si L GYang W XYang X 2009 JOSA B 26 478
8Wu YYang X 2007 Phys. Rev. 76 054425
9Wang ZShui TYu B 2014 Opt. Commun. 315 263
10Si L GYang W XLou X YHao X YYang X 2010 Phys. Rev. 82 013836
11Wang ZYu B 2014 Laser Phys. Lett. 11 035201
12Li L HLou X YLuo J MHuang Q J 2006 Phys. Rev. 74 035801
13Wang ZXu M 2009 Opt. Commun. 282 1574
14Lu X YLi J HBing L JMing L J 2006 J. Phys. 39 5161
15Wang ZChen A XBai YYang W XLee R K 2012 JOSA B 29 2891
16Xiao Y FSahin KGaddam VHua CImoto NYang L2008Opt. Express1626
17Chuang I LYamamoto Y 1996 Phys. Rev. Lett. 76 4281
18Howell JYeazell J 2000 Phys. Rev. 62 032311
19Tikhonenko VChristou JDavies B1998Phys. Rev. Lett.762698
20Doai VKhoa LBang D XHuy N 2015 Physica Scripta 90 045502
21Tian S CWan R GRong EWu J MWang HShu L JTong S LNing C ZQiang Y 2015 Physica 69 349
22Peng Yandong YJiang AMeng LLiu L 2014 Eur. Phys. J. 68 152
23Wu J LLi J B HZheng X YAnshou 2012 Opt. Commun. 285 1424
24Yavuz DSikes D E 2010 Phys. Rev. 81 035804
25Yan XWang LYin BZheng WSong HZhang JPeng Y 2008 Phys. Lett. 372 6456
26Niu YGong S 2006 Phys. Rev. 73 053811
27Niu YGong SLi RXu ZLiang X 2005 Opt. Lett. 30 3371
28Asadpour S HSahrai MSoltani AHamedi H R 2012 Phys. Lett. 376 147
29Hamedi H RJuzeliūnas G 2015 Phys. Rev. 91 053823
30Asadpour S HHamedi H RSahrai M 2012 J. Lumin. 132 2188
31Asadpour S HHamedi H REslami-Majd ASahrai M 2011 Physica 44 464
32Asadpour S HSahrai MSadighi-Bonabi RSoltani AMahrami H 2011 Physica 43 1759
33Asadpour S HSoleimani H R2015EPJ Plus301
34Asadpour S HSoleimani H R 2014 Physica 449 77
35Feili SHamedi H R 2014 Opt. Commun. 315 116
36Sadowski M LMartinez GPotemski MBerger CHeer W 2006 Phys. Rev. Lett. 97 266405
37Abergel D S LFal’ko V I 2007 Phys. Rev. 75 155430
38Ho Y HChiu Y HLin D HChang C PLin M F 2010 ACS Nano 4 1465
39Booshehri L GMielke C HRickel D GCrooker S AZhang QRen LHaroz E HRustagi AStanton C JJin ZSun ZYan ZTour J MKono J 2012 Phys. Rev. 85 205407
40Nair RBlake PGrigorenko A NNovoselov K NBooth T JStauber TPeres N M RGeim A K 2008 Science 320 1308
41Castro N A HGuinea FPeres N M RNovoselov K SGeim A K 2009 Rev. Mod. Phys. 81 109
42Yankowitz MXue J MCormode DSanchez-Yamagishi J DWatanabe KTaniguchi TJarillo-Herrero PJacquod PLeRoy B J 2012 Nat. Phys. 8 382
43Ding CYu RLi J HHao XWu Y 2014 J. Appl. Phys. 115 234301
44Ding CYu RLi J HHao XWu Y 2014 Phys. Rev. 90 043819
45Yao X HBelyanin A 2012 Phys. Rev. Lett. 108 255503
46Tokman MYao X HBelyanin A2014Phys. Rev. Lett.110077404
47Yao X HBelyanin A 2013 J. Phys.: Conden. 25 054203
48Yao X HTokman MBelyanin A2014Proceedings of SPIE 8980, Physics and Simulation of Optoelectronic Devices XXII89801A
49Asadpour S HSoleimani H R 2016 Laser Phys. Lett. 13 015204
50Hamedi H RAsadpour S H 2015 J. Appl. Phys. 117 183101
51Asadpour S HHamedi H RSoleimani H R 2015 Laser Phys. Lett. 12 045202
52Brazhnikov D VTaǐchenachev A VTumaǐkin A MYudin V IZibrov S ADudin YVasil’ev V VVelichansky V L 2006 JETP Lett. 83 64
53Jiang ZHenriksen E ATung L CWang Y JSchwartz M EHan M YKim PStormer H L2007Phys. Rev. Lett.197403
54Nesterov M LBravo-Abad JNikitin A YGarcia-Vidal F JMartin-Moreno L 2013 Laser Photon. Rev. 7 L7
55Dong HConti CMarini ABiancalana F 2013 J. Phys. B: At. Mol. Opt. Phys. 46 155401
56Cheh JZhao H2011arXiv:1107.3696
57Lukin M DYelin S FFleischhauer MScully M O 1999 Phys. Rev. 60 3225
58Scully M OZubairy M S1997Quantum OpticsCambridge, EnglandCambridge University Press
59Guo G YEbert H 1995 Phys. Rev. 51 12633
60Kosevich Y A 1997 Solid State Commun. 104 321